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Chapter 1 
Introduction 



The study of chiral symmetry breaking in quantum field theory requires the 
use of non-perturbative techniques. A commonly used method is based on 
solving the Dyson-Schwinger equations in the so-called rainbow or ladder 
approximation, which uses a bare vertex [jll], 0,0,0. However, such a trun- 
cation of the infinite set of integral equations in QED does not in general 
respect two fundamental requirements: the Ward-Takahashi identity and 
multiplicative renormalizability. This shortcoming calls in question all the 
results of the ladder approximation. 

Recently Curtis and Pennington^] presented a very interesting Ansatz for 
the full three-point vertex that satisfies the Ward-Takahashi identity and 
which is multiplicatively renormalizable to all orders in leading and next- 
to-leading logarithms in perturbation theory. They showed how powerful 
the constraints on the form of the vertex are. The natural question is to 
what extent the Curtis-Pennington modification to the ladder approxima- 
tion alters the physical results, in particular the value of the critical coupling 
for onset of chiral symmetry breaking. Applying bifurcation theory to the 
Curtis-Pennington equations, we find that there is indeed a non-zero criti- 
cal coupling; moreover the numerical changes as compared with the ladder 
approximation will be shown to be fairly small. 

These relatively minor numerical modifications encourage the hope that the 
essential results of the ladder approximation survive vertex modifications. 
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However, caution is necessary in espousing the point of view that the full 
theory has the same properties, since it is known that vacuum polariza- 
tion effects do modify the nature of the chiral-symmetry phase transition. 
Indeed, mean field behaviour appears to supervene. Nevertheless, the mu- 
tual influence of two essentially distinct but related effects, that of charge 
screening and of 'fall into the centre', have not yet been consistently taken 
into account. The basic question: 'Does the ladder approximation already 
contain the essence of the physics of chiral symmetry breaking?' is therefore 
still open. 

In this paper we shall neglect all vacuum polarization effects and we re- 
strict attention to the Landau gauge. Moreover we neglect another known 
deficiency of the ladder approximation, namely the occurrence of spurious 
first-sheet complex singularities in the fermion propagator]^. We write the 
equations from the beginning in Euclidean space and postpone the question 
of their relevance to Lorentz space. 
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Chapter 2 

Curtis-Pennington Equations 



In Euclidean space our notation for the unrenormalized fermion propagator 
isQ 

SFip)= (2.1) 
The Landau gauge bare photon propagator has the form 

D'/ik) = ^ - , (2.2) 

and the Curtis-Pennington equations are 

^(^) = i-f^ r ^\ J iy,^) (2.3) 

Art: Jo y + m'^yy) 

and 

4(.M.) = ^f -I^Jfe..), (2.4) 

The functions / and J, which depend exphcitly on m and A, were given in 
Ref. If we set m{x) = 0, we find that / and J both vanish, in fact / is 
of order while J is of order m. Hence Eq.( p.3[ ) and Eq. ( |2.4| ) possess the 
trivial solution A{x) = 1 and m{x) = 0. This is the only solution for small 
values of the coupling, oq. 

^Curtis and Pennington's S is our m and their F is our A~^. 
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In order to examine the possibility that a nontrivial solution, m{x), branches 
away from the trivial one at a critical coupling, ac, we examine the so-called 
bifurcation equation |Q. This involves differentiating Eq.(|2.3|) and Eq.(|2.4]) 
functionally with respect to m and then setting m = 0. Since / is second- 
order in m^, we find to first order in 6m that I{y,x) = and, from the 
Curtis-Pennington expressions. 



J{y,x) = 36m{y) 

3x Sm{y) — 6m{x 



y 

—9(x — y) + 0(y — x) 

X 

2 



y 



X 



y 



x^ 



6{x — y) + 6{y — a;) 



(2.5) 



According to the standard rules of bifurcation theory, we should substitute 
y for the denominator y + m'^{y) in Eq.( p.4|) . However, this introduces a 
spurious infra-red singularity, which can be eliminated by replacing the de- 
nominator hj y + m?, where the constant m can be fixed a posteriori by the 
normalization requirement 

5m{Q)=m. (2.6) 



We shall call m the dynamical mass. Rescaling x/rn^ 



y/m? 



and 6m{x)/m — > f{x), we obtain the basic bifurcation equation 



y, 



fix) 



X 

X 



ydy 
y + l 



m - f 

2 y — X 



dy j^y-^ _ x fjy) - fix) 
y + l 2 y — X 



(2.7) 



where 



A 



Sao 
47r 



In the Landau gauge, the terms involving '^^^yZi^^^ ^ire all that is left of the 
complicated Curtis-Pennington Ansatz in the linearized form appropriate to 
a bifurcation analysis. 
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Chapter 3 
Critical Coupling 



In this section we shall consider three approximate methods for handling the 
bifurcation equation (|2.7|): 

I fiy)-fix) ^ g 

y-x 

2. i^Ei^ =^ f'{y) 



fiy)- 




y- 


-X 


f{y)- 




y- 


-X 


m- 





3.1 Second-order Equation 

Here we simply ignore the extra Curtis-Pennington term. As is well-known, 
the integral equation 

n^) = -r^ny) + ^f^^f(y), (3.1) 

X Jo y + 1 Jx y + I 

then reduces to the standard hypergeometric differential equation 

xix + l)f" + 2{x+l)f' + Xf = 0, (3.2) 
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with infra-red and ultra-violet boundary conditions. The required solution is 



f{x)=2Fi{ai,a2;j;-x), (3.3) 

where 



±^\c-\; 7 = 2; Ae=i. (3.4) 

This solution of the differential equation is the only one that is finite at 
a; = 0: the normalization is fixed by the infra-red condition /(O) = 1, which 
corresponds to Eq.( |2.6| ). From the integral equation we still have to impose 
the ultra-violet condition, Bi^) = 0, where 

B{x) = ^ [xf{x)] = 2Fiiai, 0^2; 1; -x) . (3.5) 

For A < Ac, «! and 02 are real, and under these conditions the hypergeometric 
function Eq.( |3.5| ) has no finite zeros: the ultra-violet condition cannot be 
satisfied for any m 7^ 0, so there is no nontrivial solution of the integral 
equation. The imposition of the ultra-violet cut-off. A, is essential for this 
result PH , PI , 1^ , HlOl . From the physical point of view, it is easy to satisfy the 



axial- vector current conservation with a finite cut-off by just setting the bare 
mass equal to zero. Without an ultra-violet cut-off, one needs a more subtle 
treatment, distinguishing between an irregular solution, which exists for any 
coupling strength, but which leads to non-conservation of the axial-vector 
current, and the regular solution, which exists only if the coupling is large 
enough. It is not easy in a computer calculation to see the difference between 
these two types of solution unless one explicitly works with a cut-offQ. A cut- 
off is also necessary in order to make sure that the bifurcation equation is of 
classical Fredholm type, so that the standard theory can be applied^. 

When A > Ac, ai and 02 are complex conjugate numbers and the hyper- 
geometric function Eq. (|3.5|) , which remains real, has an infinite number of 
oscillations; the smallest zero, say Xq, can be chosen to define m through the 
relation 

A2 



^In the second paper of Ref.|5| Curtis and Pennington had claimed that a dynamical 
mass is generated for all values of the coupling. In fact they failed to make the above 
distinction between a regular and an irregular solution. 
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For large A , the dynamical mass is given by 



m 



A exp 



7r + ^(A) 
"2VA- Ac 



with 



(5(A) = arg 

Near the critical point, A = Ac = 

A exp 



m 



0.25, Eq.(|]6|) takes the form 
n/2 



(3.6) 



(3.7) 



v/A^ 



2 log 2 



(3.8) 



whereas the mass function (|3.3|) behaves at large momenta as 



/(x) ^logx. 



(3.9) 



The formulas (3.6) and (|3.7|) have played a great role in attracting attention 



to the non-perturbative study of QED. Firstly the critical point Oc 



has 



been interpreted 0, ^ as a nontrivial ultra-violet fixed point defining the 
continuum limit of QED. In fact, if we start from the supercritical phase 
and let A — » oo, then in order to keep the dynamical mass finite we have to 
perform the charge renormalization as follows 



a(A) 



ftc 



1 + 



log'- 

oo. Furthermore, the /5-function of the theory, 
da 2 



as A 



(3{a) = A 



dA 



a 



1 



a> an 



has an ultra-violet stable zero at a = c^c- The existence of such a critical point 
was confirmed recently by computer simulations of noncompact QED on the 
lattice Secondly, as one can see from Eg. (|3.9|) , the fermion dynamical 
mass falls off with increasing momentum as ^, and this is in sharp contrast 
to asymptotically free theories like QCD where m{jP') ~ ^ [0. It was 
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suggested ||13|| that this more slowly decreasing behaviour of the mass-function 
can solve the notorious flavour-changing neutral current problem in extended 
technicolour theories. 

By proceeding beyond the ladder approximation, we would like to find out 



whether the results ( |3.6|) and ( p.7|) are still valid or not. In Ref . [Tj] , JTS 
evidence was adduced in favour of the validity of an equation of the type 
(|3.6| ) at any finite order in the quenched approximation, and we shall make 
use of some of the arguments in Sect. 4. Our main concern here will be 
to take into account the vertex corrections in a manner compatible with 
the requirements of both the Ward-Takahashi identity and of multiplicative 
renormalizability. 



3.2 Third-order Equation 



With the second approximation, the integral equation reads 



fix) 



A ydy 



X Jo y + I 
I dy 



fiy) - If'iy) 



+A 



y + 1 



X 



fiy) - ^f'iy) 



(3.10) 



We shall show that this integral equation is equivalent to a third-order differ- 
ential equation of generalized hypergeometric type, together with an infra-red 
and two ultra-violet boundary conditions. 



A first differentiation of Eq.( p.lO ) gives 

ydy 



fix) 



JO y + I 
A /"^ dy 
2 Jx y + 1 



fiy) - p'iy) 
fiy) , 



which implies the infra-red condition 

x'fix) 



x=0 



(3.11) 



(3.12) 
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On combining Eq.( p.lOD and Eg. ( p. Ill) we find 

7 dy 



[a;/(x)]' = A 
A second differentiation gives 



[f{y) - xf\y)] . 



[a;/(x)]" 



A 



[xf\x) - /(x)] - A 



x + l 

and a third differentiation finally yields 



dy 
y + l 



f\y) 



f 



A 



X x + l 



r 



A 



x(x + 1) 



A 



x{x + 1) 



r/ = 0. 



The first ultra-violet condition (from Eq.( p.l; 



{xf{x)}' 



. a2 







(3.13) 



(3.14) 



(3.15) 



(3.16) 



is as in the second-order case; but Eg. ( |3.14| ) gives rise to a new ultra-violet 
condition, 

|[x/(x)]" + A-^(")-"^'(^)' 



x + l 



0. 



. A2 



(3.17) 



The above differential equation belongs to the class of equations of Fuchsian 
type with three regular singularities at the points x = — l,0,oo. Moreover, 
changing the variable to z = x + 1 and writing Eq. (|3.15| ) as 



z\l - z)r + [-A - (3 - A)^] zf + A/' + A/ = , 
we recognize it as the equation for the hypergeometric function, 3F2. 
We choose as a linearly independent set of solutions 



(3.18) 



m = n 



T{ai + -f^) 
r(ai - a^ + l) 



/oi +71 , Oi +72 , «i +73 
3-^2 I ai-oi+l,Oi-02+l,«i-a3+l 



(3.19) 



with the convention that the term ttj — + 1 in 3F2 is omitted when i = j. 
Here the are the roots of the characteristic equation 



+ Aa^ + (A - l)a + A = 



(3.20) 
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and 7i = , 72 = 1 + A , 73 = 2 . The a's (which should not be confused with 
the couphng, aoO evidently satisfy 



«! + ^2 + «3 = —A , aia2 + a2«3 + «3ai = A — 1 , aia2a3 = —A . 
The r factors are included for later convenience. 



The general solution of Eq. ( p.l8| ) is given by an arbitrary linear superposition 



of the fi defined in Eq.( p.l9| ). Each of these functions behaves like (1 — z)~^ 
as z — s> 1. However, N0rlund|T^ showed that the difference between any 
two of the solutions Eq. (|3.19| ) is regular at z = 1. We require this regularity 
because z = 1 corresponds to x = 0, and the infra-red condition imposes 
finiteness of /(I). We choose as two regular solutions the following: 

fiA^) = Hz)-h{z) 

/2,3W = h{z)-Uz). 



These solutions, together with any one of Eq. (p.l9|) , also form a fundamental 



system. Thus the most general solution can be written in the form 

f{.z) = 01/1,3(2;) + 02/2,3(2) + csfsiz) . 

The infra-red condition immediately gives C3 = 0, and the ultra-violet condi- 
tions lead to a system of linear homogeneous equations for the constants ci 
and C2. In order for there to be a nontrivial solution of these homogeneous 
equations, the corresponding determinant must vanish, and one obtains in 
this way the ratio cxjc^- The undetermined common constant is fixed by 
the normalization /(I) = 1. This vanishing of the determinant yields an 
equation for the dynamical mass; considering the large A limit, we find, after 
some algebra, 

i V ^ u A r(a2 + 7.)r(^3 + 7.) ^A^V"--^ 

("2 - a3){a2 + 03 - l)ai [[ ^yr? rT\ ~ 

t=i r(a2 - ttj + l)r(a3 -ai + l) J 

+two cyclic permutations of {1, 2, 3} = . (3.21) 

In the case that the a's are real, and with the ordering < 02 < ai, the 
term written explicitly in Eq.( p.21| ) is the leading one, and in this case there 
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is no nontrivial solution of Eq.( p.2lD , and hence none of Eg. ( p.7| ). However, 
on increasing A (or equivalently the couphng cto), two of the roots (cti and 
with our ordering) become equal, and then complex conjugate. For a cubic 
equation, it is easy to ascertain when this happens by looking at Cardano's 
formulas: 



as = A + 5 - ^ 



2 2 ^ ^3 



where 



A 



+ \Q 



B 



q_ 

2 



Q 



+ 



When Q < 0, Eq.( |3.2CI|) has three different real roots; but if Q > there are 
one real root and two complex conjugate ones. In the latter case two terms 
in Eq. (|3.21 ) are of the same order, and then the equation for the dynamical 
mass can be written in the form 



a2(a3 - ai)(a3 + "i - 1) H 



i=l 



r(«i + 7^) /A^ 

r(ai — + 1) \m?- 



-a 



i(a3 - a2)(a3 + "2 - 1) H 



T{a2 + ii) [A 



\ T{a2 -ai + l)\ 



^ =0. 



Writing ai = u + iv and a2 = u — iv, we obtain 

A2 



sm 

where the phase is given by 

(5(A) = ar 

and the mass is 

m = A exp 



flog ■ 



5(A) 



a2{a3 - ai){a3 + ai - 1) [[ 

i= 



T{ai -ai + 1] 



Tin + S{\) 



2v 



n = l,2. 



12 



As usual, the value n = 1 corresponds to the ground state of the system. 
The critical coupling constant Ac(= separating the chirally symmetric 
phase from that of spontaneous chiral symmetry breaking, is determined by 
the characteristic equation Q — 0, i.e. 



(3.22) 



3A^ - 12X^ + 32A2 + 12A - 4 = . 
Numerical calculations give Ac = 0.2172 and 

ai(Ac) = a2{Xc) =Uc = 0.4435 , = -Ae - 2uc = -1.1042 . 

The above value of Ac implies Oc = 0.9098. 

Near the critical value A ~ Ac, the imaginary part of the roots behaves like 



v{X) ~ \/A - Ac , 

while S{X) ~ V. Finally we obtain an expression for the dynamical mass: 



m = A exp 
where the constants (Xj are given by 

71 p{Xc) 



+ (72 



(3.23) 



(7i = 
(72 = 



2 [3g'(Ac)]^ 



1 1 

— + ^ 27 

Uc 1 + Ac + Mc 

+ '0(Ac + 3uc) - i>{l + Ac + Uc) - ip{2 + Uc) - ip{uc) 



(3.24) 
(3.25) 



Here ip{x) = dlogT{x)/dx, 7 is Euler's constant. Numerically, we find 

(Ji = 1.525 (J2 = 1.604 . (3.26) 
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3.3 Fourth-order Equation 



Our third approximation to the Curtis-Pennington contribution leads to the 
following fourth-order differential equation: 



A 



2 

where 



1 + A 



X 



x + 1 



Xf + 



A x^ 



2 (x + 1)2 



(3.27) 



(j){x) 



y + 1 



jx'^ — X + log{x + 1) ~ 



(3.28) 



Substituting / ~ x we obtain the characteristic equation 

- 2(2 - \)a^ - SXa"^ + 4(1 - A)a - 4A = , 

with the critical value Ac = 0.23320 and the roots 

"1 = "2 = 0.4851 

as = -1.1113 
a4 = 15.294. 

These a's and also ai can be calculated using the explicit form for the roots 
of Eq.(^). 
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Chapter 4 
Discussion 



A comparison of the results in the third- and fourth-order cases with those in 
the second-order case [Eg. ( p. 81 )] shows that the numerical changes are rather 
small (recall that = ^)|]: 







Uc 


0"! 




2nd order 


1.047 


0.5 


1.571 


1.386 


3rd order 


0.910 


0.443 


1.525 


1.604 


4th order 


0.980 


0.485 


1.462 





We propose to show now that the form Eq.( |3.23| ) is generic; that is, at any 
order of the quenched approximation, the same asymptotic formula is valid, 
with suitable values of the parameters Ac, and a2- Holdom's |14] essen- 
tially diagrammatic analysis, within the quenched approximation, leads to a 
bifurcation equation for the mass function of the following form 

fix)=\ f^dyy^K[x,y;0], (4.1) 
Jo 1/ + 1 

where the kernel K[x,y;0] is symmetric under x ^ y and xK[x,y]0] is a 
function of x/y only. This is known to be in agreement with a perturbative 

^ These numbers for ac are in agreement with very recent numerical resuhs of Curtis 
and Pennington 
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expansion to order 0(a^)[|T8|, and it leads immediately to a dynamical mass 



of the form Eq.( 3.23 ). 



Here we have shown that the same result also holds when we use a non- 
perturbative Ansatz for the vertex based on the Ward identity and multi- 
plicative renormalizability, although this vertex does not lead to a bifurcation 



equation of the form Eq. (|4.1| ). In general this result will hold for all approx- 
imation schemes in which the integral equation for the mass can be reduced 
to a differential equation of the Fuchsian type. 

Once we have a differential equation, possibly of very high order, it is easy to 
understand that the behaviour of the mass function at large momenta is 
where a is a root of an (algebraic) characteristic equation of finite order (say 
n), with real coefficients. The onset of oscillations corresponds to equality 
of two of the roots of this characteristic equation: for larger values of the 
coupling, the two roots are complex conjugate and oscillatory behaviour is 
manifested and the ultra-violet conditions can be satisfied. To calculate the 
imaginary part v of these roots, let us consider the characteristic equation's 
discriminant: 

A=«o""' n K 

l<i<j<n 

(ao is the coefficient of the largest power, a"). The discriminant is symmet- 
ric with respect to the roots a^, and it may be expressed in terms of the 
coefficients of the characteristic equation (for the rules to calculate the dis- 
criminant, see for example |T^). Thus A is an explicit function of A. Note 



that A changes sign when a pair of complex conjugate roots appears. 
The critical couphng constant Ac is determined by the condition 

A(Ae) = (4.2) 

and the imaginary part v vanishes like A^. Normally Ac is a simple root of 
Eq.(|]|), so that 

^_A'(Ac)v/A-Ac, 

and this leads to formula Eq.( p.23D for the dynamical mass. In principle Ac 
could be a degenerate root, but this would be an 'accident': the generic case 
is as just given. 
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Our considerations show that the dependence of the dynamical mass on the 
couphng constant, Eq. (|3.23|) , is a universal feature of quenched theories. The 
behaviour of the real part of those roots, u, is however not universal. It 
gives the asymptotic power-envelope of the decreasing and oscillating mass 
function: 

2-7m 

where 7^ is the anomalous dimension of the composite operator ipip. In the 
lowest approximation (the second-order equation), we have 7^ = 1, but the 
Curtis-Pennington Ansatz leads to a modification of this dimension. For 
example, in the cubic approximation we find 7^ = 1.114. Values 7m > 1 
imply that some four-fermion operators like {ip'ipy, become relevant 

in the continuum limit and must be added to the QED Lagrangian^], which 
means that pure QED is not self-consistent. 



1 


u 
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